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ABSTRACT  We  have  developed  a  crystallographic-based  viscoplastic  model  for  non- 
isothermal  high  temperature  deformation  and  coupled  it  with  damage  kinetics.  Several 
damage  mechanisms,  namely  the  multiplication  of  mobile  dislocations,  void  growth  and 
including  scale  effects,  caused  by  dislocation  extrusions/intrusions,  have  been  considered. 
We  applied  two  body  interaction  concepts  from  chemical  kinetics  to  include  in  our 
constitutive  relations  the  generation  and  interaction  of  pinned  and  mobile  dislocations. 
The  relative  importance  of  each  of  the  damage  modes  at  different  temperatures  and  stress 
levels  has  been  deduced  by  comparing  strain-stress  and  crystallographic  texture 
predictions  calibration  against  test  data  and  has  shown  that  the  dislocation  density  growth 
may  be  important  in  characterizing  primary  creep. 

INTRODUCTION:  An  understanding  of  the  micromechanics  of  high  temperature  creep 
and  damage  accumulation  in  single  crystal  nickel  base  superalloys  is  important  for  the 
design  of  turbine  blades  and  vanes  in  advanced  commercial  and  military  gas  turbines. 
Our  models  must  not  only  predict  the  creep  response  of  materials  in  gas  turbines  they 
must  also  be  able  to  model  virgin  material  since  most  experimental  data  has  been 
obtained  from  “as  received”  materials  and  “as  received”  materials  have  considerably 
smaller  dislocation  densities  than  pre-strained/pre-exposed  materials.  This  can  drastically 
change  the  primary  creep  response.  The  Orowan  equation,  where  strain  rate  is 
proportional  to  the  dislocation  density,  is  commonly  used  to  describe  the  mechanical 
response  of  materials.  But  clearly  if  the  dislocation  density  increases  then  the  strain  will 
also  increase  and  the  equations  must  include  an  additional  tenn  proportional  to  rate  of 
growth  of  dislocation  density.  We  believe  that  including  dislocation  density  growth  rate 
in  our  model  may  predict  the  observed  experimental  differences  between  “as  received” 
materials  and  pre-strained  materials. 

PROCEDURES,  RESULTS,  AND  DISCUSSION:  We  have  used  a  standard 
viscoplastic  power  law  creep  with  a  back  stress  to  represent  the  response  of  the  material. 
The  constitutive  law  for  the  inelastic  strain  rate,  y?  along  slip  plane  i  is  proportional  to 
the  mobile  dislocation  density,  pm ,  and  is  written  as  (Staroslesky  et  al,  2006) 

yP  =  y0(pm  / p0) / s*\  sgn((ri  -®i)/5*j ,  where  y0  is  a  time  constant,  xi  is  the  slip 

plane  resolved  shear  stress,  s*  is  the  isotropic  yield  stress,  0)l  is  the  slip  plane  back  stress, 

•  •  •  ^  • 
and(  ')is  the  rate  of  change  with  respect  to  time.  Both  st  and  a>i  have  evolution 
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equations  of  their  own.  We  obtain  relations  for  mobile  and  pinned  dislocations  along  each 
slip  system  as: 


Pm=  M 
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where  M  and  n  represents  specific  time  constants,  different  for  the  octahedral  and  cube 
slip  systems,  p “  ,  is  the  saturated  mobile  dislocation  density,  psps ,  is  the  saturated  pinned 


dislocation  density,  and  s2  is  a  positive  constant,  vis  the  pore  volume  fraction  and  has  an 
evolution  equation  associated  with  it.  Please  see  Staroselsky,  et  al  (2006)  for  a  more 
complete  discussion  of  the  model.  The  Orowan  (1940)  equation  is  commonly  used  to 
describe  the  mechanical  response  of  materials  and  is  referenced  widely  [e.g.  see 
Thompson,  2006], 


The  Orowan  equation  is  generally  written  for  edge  dislocations  with  a  dislocation 
density,  p,  as  y  =  pbv ,  where  y  is  the  shear  strain,  v  is  the  dislocation  velocity,  b  is 
the  Burgers’  vector,  and  (’)  =  d(  )/ dt .  However  if  the  dislocation  density  increases  then 
the  strain  will  also  increase  and  the  equation  should  include  a  p .  For  example,  we  could 

write  y  =  pbL  +  pbL  .  The  quantity  L  is  some  characteristic  dimension  which  will  make 
the  strain  rate  the  total  derivative  of  a  function.  An  equation  containing  p  is  not  a  new 
idea,  and  has  been  proposed  independently  by  others  [e.g.,  see  Gupta,  et  al,  1974].  The 
characteristic  distance  L  is  the  key  parameter  that  must  be  developed.  Our  proposal  is  to 
assign  to  L  the  distance  between  pinned  dislocations  similar  to  Thompson  (2006)  where 
the  author  has  proposed  using  the  average  dislocation  displacement. 


For  example  consider  the  singular  components  of  the  strain  for  edge  and  screw 
dislocations  describe  the  state  of  strain  in  a  material.  For  an  edge  dislocation  in  the  x-y 
plane  (i.e.,  the  dislocation  line  is  parallel  to  the  z-axis),  and  the  discontinuous  slip  along 
the  x-axis,  the  singular  part  of  the  displacements  in  cylindrical  coordinates  (Nabarro, 
1987)  are,  u  =  b/2;r{#,ln(Wb),0} ,  where  the  vector  components  are  along  the 
r,6,z  are  the  directions  respectively.  For  a  screw  dislocation  in  the  x-y  plane  the  singular 
part  of  the  cylindrical  displacements  (Nabarro,  1987)  are  u  =  b!2n{ 0,0,6*} .  For  the  edge 

dislocation  the  small  strain  components  can  be  evaluated,  and  then  the  maximum  shear 
strain  can  be  readily  shown  to  be  <?max  =  b  /  2nr  in  cylindrical  coordinates.  Similarly  for 

a  screw  dislocation  the  maximum  shear  strain  is  sm„  =  bl  Anr .  Recall  that  the  maximum 

IIldA 

engineering  shear  strain,  ymax  =  2smax , 
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For  dislocations  that  are  distributed  throughout  an  area  with  density,  p ,  the  strain  tensor 

becomes  for  both  edge  and  screw  dislocations  can  be  written  as  s y  =  \pbfij(r~R)dS , 

s 

where  S  is  the  area  over  which  the  dislocations  are  distributed,  and  we  have  assumed 
that  the  rate  of  change  of  the  integrand  over  time,  t,  is  slow  enough  for  a  steady  state  in 
the  strains  to  exist.  The  dislocation  density  can  include  both  edge,  pe ,  and  screw 


dislocation  densities,  ps  ,  where  p  =  pe  +  ~^Ps  •  The  strain  rate  can  be  found  by 
differentiating  with  respect  to  time,  then  Gy  =  J  \pbfij  -  pbR  •  V fij  \lS  ,  where 
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(')  =  3(  )/dt,  f  =  f(r-R ),  and  r  is  the  vector  (x,y,0)  while  R  is  the  vector 
( X ,  7,0) .  The  integration  is  perfonned  over  the  coordinates  (x,  y,0) .  If  p ,  p ,  and  R 

are  constant  throughout  the  area  S  then  s-  =  />/)  j  fplS  -  pbR  •  j  fpidC  where  the 

s  c 

contour  C  surrounds  the  area  S  ,  and  use  has  been  made  of  the  divergence  theorem. 


We  can  find  an  approximate  value  for  the  second  integral  by  considering  the  case  of  a 
circular  contour  of  radius  a  surrounding  the  area  S  ,  and  all  of  the  dislocations  moving 
with  velocity  v .  Then  the  maximum  engineering  strain  rate  for  edge  dislocations  with 
p  =  0  is  given  by  y|Tiax  =  pebv ,  which  exactly  matches  the  Orowan  equation.  We  now 
need  an  estimate  for  the  term  proportional  to  p .  Consider  again  the  case  of  a  circular 
area  of  radius  a  .  An  integration  yields  y|Tiax  =  2 peba  +  pebv .  Our  proposal  is  to  make 
the  tenn  2 a  in  the  p  term  proportional  to  the  average  distance  between  pinned 
dislocations.  That  is,  if  the  pinned  dislocation  density  is  p  then  the  distance  is 

proportional  to  1/  ^Jp^ . 
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